Minimal model of self-replicating nanocells: a physically embodied 
information-free scenario 
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The building of minimal self-reproducing systems with a physical embodiment (generically called 
protocells) is a great challenge, with implications for both theory and applied sciences. Although 
the classical view of a living protocell assumes that it includes information-carrying molecules as 
an essential ingredient, a dividing cell-like structure can be built from a metabolism-container 
coupled system, only. An example of such a system, modeled with dissipative particle dynamics, 
is presented here. This article demonstrates how a simple coupling between a precursor molecule 
and surfactant molecules forming micelles can experience a growth-division cycle in a predictable 
manner, and analyzes the influence of crucial parameters on this replication cycle. Implications 
of these results for origins of cellular life and living technology are outlined. 
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I. INTRODUCTION 

The transition from non-living to living systems cov- 
ers a b road spectrum of increasi ngly complex organi- 
zation l(Smith fc Szathmarvl llQQfj) . One of such first 
steps in this transition leads from ordinary chemical 
auto-catalysis to self-replication. The latter can be dis- 
tinguished from the former by the existence of self- 
bounded entities which produce copies of themselves, 
rather than a mere increase in chemical concentration. 
Therefore, self-replication relies on organization prin- 
ciples unlikely to be found in homogeneous solutions. 
Within the last years increasing attention has been payed 
to the possibility of building small-scal e protocells, in 
parti cular using a bottom-up approach l|Szostack et al.l 
l200lj) where the building blocks (not necessarily from 
biotic origin) would assemble spontaneously and, un- 
der appropriate conditions, develop a growth-fission cy- 
cle. While extensive research has been perfo rmed on the 
self-reproducing capabilitie s of bio-polymers ijKiedrowskil 
1986, iTiivikua et al"] ll990b even much simpler systems 
can be driven into dynamics that we identify as self- 
replication. In this context, micelles have been pro- 
posed to serve as life-like structures a ble to undergo self- 
reproduction ijBachmann et al.1ll992j) . 

These micellar systems can be considered to be sim- 
pler in organization than bio-polymers because they lack 
any genetic information that could be passed from one 
generation to the next. In this context, protocellular en- 
tities lacking information woul d be under th e umbrella 
of Oparin's views of life origins l)Oriarinl [l936). who sug- 
gested that primitive self-replicating vesicles would have 
predated information-based cells. Similarly, other au- 
thors have advocated for this scenario under a computa- 




FIG. 1 The basic model of nanocell replication explored 
in this paper. Here small-sized micelles are formed by am- 
phiphiles (here indicated as H — T connected pairs of balls). 
These amphiphiles have a hydrophilic head (H) and a hy- 
drophobic tail (T). Precursor molecules are also shown as 
two connected, smaller open balls, both of them of hydropho- 
bic character. Under the presence of catalyically active am- 
phiphiles, precursors are transformed into additional am- 
phiphiles. Incorporation of many such building blocks allows 
the nanocell to grow in size. When a critical size has been 
reached, the nanocell looses its stability and divides into two 
smaller aggregates thereby closing the replication cycle. 
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tional perspective. Under this view, Dyson ljDvsonlll999j) 
indicated that current cells involved both software (DNA 
and RNA) and hardware (protein machinery). Although 
hardware can exist under the absence of software, the op- 
posite is not allowed to occur. We can add to this picture 
of cells that the presence of a container is a very impor- 
tant piece for the hardware to properly work: only when 
the pieces are close together we can expect the machinery 
to operate. 

The nanocellular system discussed here is based on a 
micelle that is coupled to a minimal metabolism (figure 
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The system is constantly driven away from equilib- 
rium by the supply of precursors, which are supposed 
to have higher internal energy than the surfactants that 
build up the micelles. The metabolism transforms pre- 
cursor molecules into new building blocks of the nanocell. 
The metabolic turnover is thereby enhanced by its own 
outcome or — in other words — the metabolism is an auto- 
catalytic turnover of precursors into new surfactants. 
This system rese mbles one studied exp erimentally by 
Bachmann et al. (Bac hmann et al.l ll992). However, in 
their work the catalytic effect is a feature of the micelles 
(probably due to milieu effects), while, in our system, it 
is a feature of the molecules themselves. 

Either way, the replication mechanism of micelles en- 
visioned both in Bachmann's as well as our system is 
the following: micelles incorporate hydrophobic precur- 
sor molecules where they are afterwards transformed into 
new surfactants. Due to this process, the number of sur- 
factants increases, while at the same time, the volume 
of the hydrophobic core becomes smaller. It is assumed, 
that when a critical ratio of surfactants versus core vol- 
ume is passed, the aggregates become unstable and will 
divide into two daughter cells. Whether or not the exper- 
imental system follows this pathway has, to our knowl- 
edge, never be clarified. However, once such a replication 
mechanism is considered the basis for a more complicated 
protocell design, knowledge of the exact replication kinet- 
ics becomes crucial. 

Together with experimental approaches aimed to find- 
ing the conditions for protocell replication to occur, there 
is an increasing need of appropriate, well-grounded theo- 
retical and computational models. Such models can help 
understanding the constraints that might operate in the 
self-assembly of micelles and other molecules and how 
they can properly trigger growth and splitting. Our work 
is a first step in this direction. 



II. DPD NANOCELL MODEL 

In this paper we consider a dissipative particle dynam- 
ics (DPD) approach to modeling embodied protocells 
employing a physically and chemically simplified formal- 
ism. This is one possible approach that can be used 
in order to tackle the complexities of molecular aggre- 
gates. Other approaches include Molecular Dynamics 
(MD), Brownian Dynamics (BD) and Monte Carlo al- 
gorithms. Each o f these methods has i t s own advantages 
and drawbacks llBinder fc Heermannl 119971 IHeermannI 
11990^ and a compared analys is has been presented else- 
where llChen fc Szostaldl2004ft . 

Previous work done by Ono a nd Ikegami involved a 
lattice-based, protocell dynamics ijOno fc Ikegamilll999|) . 
These structures display a special type of cell-like replica- 
tion. While they are remarkable in their self-organizing 
behavior, the underlying rules of interaction and the spe- 
cial properties of the membrane-like structures restrict 
their relevance to the arena of artificial life. By using 



more appropriate molecular interactions within a three- 
dimensional, water-filled environment (see for example 
Rasmussen an d Nilsson's cellular auto mata approach to 
self-assembly llRasmussen et al1l2003|) l we seek to pro- 
vide the (as far as we know) first computational evidence 
that such a simple protocellular cycle is feasible. 

A. Dissipative particle dynamics 

DPD is a coarse grained, particle based simula- 
tion technique comparable to Brownian Dynamics. 
It was proposed by Hoogerb rugge and Koelmann 
ijHoogerbrugge fe KoelmarJll992() an d gained significant 
theoretical support in the late 1990s jEs panol fc Warrenl 
119951 iGroot fe Warrenlll997llMarshlll998(l . In the con- 
tcxt of biological systems, DPD models have been suc- 
cessfully used to captu re the dyn amics of membranes 
llVenturoli fc Smitlll999l). ve sicles l)Yamamoto fc Hvodd 
| 2003t lYamamotq et alJl2002ft and micelles l|Groodl2000l 
lYuan et aUl2002|) . 

A DPD simulation consists of a set of N particles which 
are described by their type, mass m,-, position r,, and 
momentum = TOiVj. These particles — usually called 
beads throughout the literature — are not meant to repre- 
sent individual atoms. Instead, they represent groups of 
atoms within a molecule (like several CH2 groups within 
a hydrocarbon chain) or even a group of small molecules 
such as water. 

Newton's Law of motion is used to determine the tra- 
jectory of each individual bead: 



The force Fj that acts on particle i is expressed as the 
superposition of pairwise interactions 

N 

F s : = 5>y ( 2 ) 

3=1 

In Newtonian dynamics, the (central) force FR can be ex- 
pressed as the negative gradient of a potential 0y , namely 

*ii = -V0y (3) 

The resulting dynamics are conservative and obey the 
Hamiltonian 

1 N 1 N 

H =2^ miVl + 2 ^ (4) 

i=l i,i=l 

While this approach is undertaken in molecular dynamics 
simulations, coarse grained simulation techniques try to 
aggregate some of the molecular degrees of freedom by 
the use of the so called Langevin formalism: additional 
forces F,^ and FfJ are added to the conservative force 
to express friction and thermal motion. They introduce 
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energy flows between the explicitly modeled mesoscale 
and an the underlying microscale. Together, they act as 
a thermostat to regulate the effective temperature, i. e. 
mean velocity, of the system. 

The thermostat used in the DPD formalism is given 
by the equation 

F S + V u = ('/-"-., "»<, • v tf ) + o-w 2 (r tf )^ tf ) n,„ (5) 



where rv. 



r,-| is the distance, n 



(n - Tj)/n. 



ij 



Vj the relative ve- 



the (unit) direction, and v 
locity between beads i and j. 77 is the friction coefficient 
of the fluid and a the amplitude of thermal motion. £y 
is a random variable with £y = , and otherwise Gaus- 
sian statistics, w is a distance weighing function usually 
defined as 



(6) 



w(r) = 21- 



where r c is a cutoff range used to limit the maximal range 
of all interactions for performance reasons. It can be 
shown that the equilibrium temperature of the syst em 
tends towards k b T = a 2 /2-q ijEspaiiol k, Warrenlll995|) . 

While there is a variety of other thermostats used in 
coarse grained particle simulations, the above mechanism 
is unique in that it both conserves linear and angular mo- 
menta and fulfills the fluctuation-dissipation theorem. As 
a consequence of the former, the resulting dynamics are 
consistent with the Navier- Stokes equations and hence 
preserve laminar flow properties of the system. The lat- 
ter property ensures an energy distribution in the system 
following Maxwell-Boltzmann statistics. The overall dy- 
namics, therefore, capture both hydrodynamic and ther- 
modynamic traits of the systems. 

In almost all DPD studies, the conservative force is 
derived from a soft-core potential of the shape 



<pij (r) 



M) ; 





if t < r c 
if r > r c 



(7) 



The potential energy expressed by <j>ij should not be un- 
derstood as the mechanical energy, i. e. enthalpy, of 
the system, but rather as its f ree energy contribution 
l|Pagonabarraga fc Frenke]| l2001'). Following this ratio- 
nale, the interaction parameters are used to express 
dissimilarities of substances due to high enthalpy as well 
as entropy contributions, respectively. Therefore, they 
can be related to Flory-Huggins coefficients known from 
polymer theory. 

For the study of lipids and surfactants, covalent bonds 
between beads are commonly introduced as harmonic 
spring forces: on top of the above interactions, bonded 
beads interact according to the potential 



I J 



(8) 



where b is the strength and r b the optimal distance of 
covalent bonds. As usual, we use r c , to, and k b T as units 



of space, mass, and energy, resp ectively. The time unit 
follows from equation ^ as r = ^Jm/k b Tr c . 

To model the system under consideration, we define 
beads of type W (water), H (hydrophilic "heads") and 
T (hydrophobic "tails" of amphiphile s) with interaction 
parameters taken from l|Grootl 1200(1) . unless otherwise 
specified: 



W 
H 
T 



W 



H 



25k b T 15k b T 80k b T 
lbk b T Zbk b T 80k b T 
80fc b T 80k b T 15k b T 



All beads have mass Ito. Precursor molecules are mod- 
eled as dimers of bonded T beads, surfactants as dimers 
of one T and one H bead. Here we have: b — l25k b T, 
r b = 0.5r c for all covalent bonds. 

The objective behind this parameter set is to model 
surfactants that form spherical micelles. To achieve this, 
the effective head area must be large compared to the 
volume of the hydrophobic core (packing parameter 1/3). 
This is expressed by ott < »ww < hhh- Furthermore, 
surfactant heads have a high affinity to water («hw < 
aww)i which is usually due to charges in the hydrophilic 
groups of the molecules. This assumption ensures that 
aggregates with high surface area (spherical micelles) are 
prefered over aggregates with less surface area (rod-like 
micelles) in the process of total energy minimization. 



B. Incorporation of chemistry 

The metabolic reaction under consideration takes the 
following form 



T-T 



H T 



(9) 



This reaction is modeled by a stochastic process that has 
formerly b een used in Brownian Dynamics simulations 
llOnoll200lD . Inbetween every two steps of the numerical 
integrator for the DPD equation of motion, each precur- 
sor dimer can be transformed into a surfactant molecule 
with a spontaneous reaction rate k b . The spontaneous 
reaction can be catalytically enhanced by nearby surfac- 
tants whose catalytic influence decreases linear with the 
distance to the reactant up to a certain threshold r cat . 
For simplicity, the effect of several catalysts is modelled 
as a superposition: 



i-SJ iir<r ^ (io) 
otherwise 



where rc is the distance of the catalyst and k s the maxi- 
mal catalytic rate per catalyst. For clarity of the results, 
we set the spontaneous reaction rate in our simulations to 
Or -1 . In the upcoming simulations, k cat is set to l.Or -1 , 
r cat to lr c . If a reaction occurs, the type of one random 
T bead is changed to H, but positions and momenta are 
preserved. 



4 



We also introduce particle exchange into the model 
to mimic the support of new precursors into the system: 
During the simulation, two water beads can be exchanged 
by precursor dimers with the probability 2.5T0~ 5 r~ 1 per 
water bead within a region of radius 2r c . Again, bead 
positions and momenta are preserved. 



III. RESULTS 

We implemented the DPD method using a velocity- 
Verlet integrator (A = 0.5) with a step width of 0.04r. 
The spatial domain is three dimensional, with periodic 
boundary conditions and size (10r c ) 3 . In all the follow- 
ing simulations, the system is initialized with one sur- 
factant dimer and 2998 water molecules adding up to a 
mean particle density of 3r~ 3 . Unless otherwise specified, 
simulations have been run for Or < t < lOOOr (25000 it- 
erations). 

For simulation runs with the parameter set given in 
the last section, we typically get the following behav- 
ior: water beads are successively exchanged by precur- 
sors in the exchange region of the system. While diffus- 
ing through the system, precursors form droplets due to 
their hydrophobic trait. Once such a droplet reaches the 
initial surfactant, the precursors are turned into surfac- 
tants. The turnover happens fast compared to supply 
and diffusion of the precursor. New surfactants quickly 
arrange into a micellar shape with hydrophobic beads in 
the interior and hydrophilic beads towards the surface 
of the assembly. With the parameters introduced be- 
fore, this rearrangement takes only few time steps and 
is thus fast compared to the metabolic turnover. Such 
spontaneously formed nanocells diffuse through the sys- 
tem space as aggregates and eventually incorporate ad- 
ditional precursor droplets in their interior, where the 
metabolic process is repeated. The evolution of the over- 
all system composition (number of beads per type) traces 
the different processes on their respective time scales (see 
figure |21 for bht = awt = 80fcf,T (upper panel) and 
ttHT = «WT = 120/cfcT (lower panel). As one can see, 
the overall production of surfactants is limited by a linear 
growth that results from the constant supply of precur- 
sors. Locally, however, when a single droplet is consumed 
by a nanocell, the metabolic turnover exhibits a logistic 
growth which is suspected from resource limited auto- 
catalysis. The logistic growth can be best seen in the 
lower panel of figure [21 around t = 700r. 

On the level of individual micelles, the metabolic pro- 
cess increases the aggregate number, i. e. the number of 
surfactants per micelle. In a pure water-surfactant sys- 
tem, micelles would reject the surplus of surfactants into 
the bulk phase. In an oil-surfactant-water system, as the 
one under consideration, the hydrophobic core formed by 
the precursors, stabilizes the assembly far beyond its orig- 
inal aggregate number. As a consequence, we could ob- 
serve that nanocells increase in aggregate number when 
new surfactants are synthesized. While the precursor sur- 
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FIG. 2 Evolution of bead numbers with time for two different 
hydrophobicity values (oht and owt). The constant supply 
of precursors is counteracted by their transformation into sur- 
factants. This transformation happens in spurts rather than 
continuously, as the precursor forms droplets in the aqueous 
solution. Size and frequency of these spurts depend on the 
hydrophobicity of the tail beads. 



factant ratio shifts, the nanocell changes its shape from 
a spherical to a rod-like micelle. When all or nearly all 
of the precursor is turned into surfactant, the nanocell 
finally becomes unstable and divides into two smaller 
aggregates (see figures |3 and 0J. The nanocell division 
occurs in the cylindrical middle-part of the rod-like ag- 
gregate by indentation of surfactant heads. It induces 
vibrating modes into the daughter aggregates as they re- 
arrange back to a spherical shape. Although this vibra- 
tion is successively dissipated into undirected motion, it 
sometimes leads back to short series of temporary fusion 
and fission of the daughter cells. 

Elongated micellar structures are well-known from 
worm-like micelles which usually consist of two surfac- 
tants with different curvature. Such worm-like micelles 
are stable equilibrium a ggregates and exhibit an e xpo- 
nential size distribution ijKroger fc M akhloufi Il99"r^ . In 
the system studied in this work, however, the elongated 
aggregate is not stabilized by different curvature values 
of the components, but by the hydrophobic core. Ac- 
cordingly, once the precursor is turned over into new 
surfactant molecules, the elongated structure looses its 
stability. 

There is a second pathway, however, that might jeopar- 
dize the above scenario: once in a while throughout our 
simulations, nanocells loose individual surfactants into 
the bulk phase. If this relaxation process happens fast 
compared to the metabolic turnover, the nanocells might 
not be able to reach the division size. Surfactants in the 
bulk phase may however metabolize precursor droplets 
and spontaneously form nanocells on their own. 

Formally, fission events can be written in the form of 
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FIG. 3 Here, the size evolution of figure [5] (lower panel) is 
tracked for individual nanocells: each line designates the size 
evolution of a single aggregate. Horizontal lines result from 
fission (or dissociation) events, after which two lines indicate 
the fate of the daughter cells. Isolated dots denote short term 
vibrations during which nanocells divide and fuse within less 
than It. For clarity, such horizontal lines have been sup- 
pressed for such vibrations. As one can see, only two daughter 
cells result from a true fission event (at t — 722r). The other 
two result from dissociation of single surfactants that start 
to turn over precursor droplets found in bulk phase. These 
surfactant dissociations happen at t — 324r and t = 684r. 
Furthermore, two nanocells fuse at t = 655r. 



as well as fissions into and fusions from more than two 
aggregates. For all the simulations performed, such out- 
liers amount to less than 0.3% of the total transitions. 
They have been neglected for further analysis. 

In the remainder of this work, we will analyze two key 
parameters of the model and their influence on the dy- 
namics of the system. Thereby, we will put our attention 
on the fission pathways discussed above. The first pa- 
rameter we analyze is the hydrophobicity, i. e. the values 
otw and ath- This parameter is related to the dis- 
similarity between surfactant tails and water. For most 
surfactants the hydrophobicity is solely a function of the 
lengt h of the hydrocarbon chains ijEvans &: Wennerstrdml 
1999). The value is thus easily adjustable in experimental 
setups. 

Second, we analyze the influence of the catalytic rate 
k c — hence, the speed of metabolic turnover — on the di- 
vision pathway of the nanocell. As we did not specify 
the molecular implementation of precursor and surfac- 
tant, it is conceivable that in an experimental setup, one 
can choose from a wide range of catalytic efficiencies. 



a chemical reaction: 



Sn+m 



~^ S n + S n 



(11) 



where S n , S m and S n + m are aggregates of size n, m, 
and n + m, respectively, and fc~ n is the fission rate. For 
m = 1, one obtains dissociations as a special case. Anal- 
ogously, association and fusion events take the form 



S n 



+m 



(12) 



In order to quantify micellar fission and surfactant dis- 
sociations, nanocells have been identified by a variant of 
the flood fill algorithm: every two T beads within a dis- 
tance of lr c or less have been considered to belong to 
the same aggregate. The aggregate number is defined as 
the number of participating surfactants. This allows for 
tracking individual aggregate sizes and their transitions 
through time. Each of the above reaction schemes results 
in two transitions, given by 



S n +r 



and 



Sn 



s 



n+m 



(13) 



(14) 



Furthermore, the turnover of precursors results in the 
transition 



Si 



ra+l 



(15) 



Not all transitions can be expressed by the chemical 
reaction scheme given by ^2 and Those transitions 
are of the type 



S 



n+m 



Si 



S n + S m +l 



A. Influence of hydrophobicity on nanocell dynamics 

For three hydrophobicities, 40, 90 and 120fcf,T, his- 
tograms of such transitions are shown in figure [S] The 
figures reveal a clear trend both in aggregate numbers 
as well as transition types. For hydrophobicity A0k{,T, 
the system is almost entirely composed of single surfac- 
tants and small aggregates in bulk phase. 33.4 % of the 
transitions arc dissociations and associations of two sin- 
gle surfactants. For the few bigger aggregates, transi- 
tions are distributed more or less homogeneously, i. e. 
surfactant dissociation is as likely as proper aggregate 
fission. Thus, for weak hydrophobicities, the system re- 
sembles a homogeneous solution without significant for- 
mation of structures. For 120fcf,T, on the other hand, the 
transition histogram looks completely different. Associ- 
ations or dissociations of two isolated surfactants make 
only 3.7% of the transitions, for this parameter. The 
most prominent transition type is the turnover of a pre- 
cursors within nanocells that range in size from 1 to 25 
surfactants (20.6% of all transitions), represented by high 
values in the lower secondary diagonal. The absence of 
an upper secondary diagonal reveals that there are no 
surfactant dissociations except some between aggregate 
numbers 8 and 16. Due to the higher stability of aggre- 
gates, there are far less overall transitions than in the 
previous case (321 compared to 24743), which poses a 
problem when trying to obtain smooth histograms. Al- 
though there are distinct cases of proper fission events 
(e. g. -* Su + Sn, S 20 -* Sn + S 9 ) it is hard 

to tell from the histograms whether such fission is more 
likely than single surfactant dissociation. Varying the 
hydrophobicity between these two extrema yields traits 
of both other histograms. One example is given in figure 
El for 90/c h T. 



FIG. 4 Metabolism and fission of a nanocell (T beads are shown in light, H beads in dark gray — water not shown): The 
initial surfactant metabolizes a precursor droplet and turns it into a functioning nanocell (panels 1-3). While the precursor is 
consumed, the nanocell elongates to account for the changing precursor surfactant ratio (panel 4-5). Such elongated structures 
can be stable for several time units, until — when all precursors are turned into surfactants — the nanocell divides into two 
daughter cells (last panel). 
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FIG. 5 Histograms of nanocell size transitions for different 
hydrophobicities (see text on how transitions are defined). 
For a transition S n — » S m , n is ordered along the vertical, m 
along the horizontal axis. Colors indicate the number of oc- 
curences. Note that colors have been scaled by root functions 
to emphasize seldom events. 



As new precursors are constantly supplied and 
nanocells grow and divide over time, it is somewhat diffi- 
cult to capture mean aggregate numbers of the assem- 
blies. Nevertheless, these values are prominent char- 
acteristics in the study of micellar systems and their 
knowledge can help to get insight into the system under 
consideration. To gain aggregate numbers we compute 
the average value m + n for all transitions of the form 
Sm+n —> S n + S m , i. e. we average aggregate numbers 
only in moments of fission or dissociation. The results 
can be regarded as the mean maximal aggregation num- 
ber of nanocells. Figure shows results for hydrophobic- 
ity values from 40fc&T to 120/cfcT together with a simple 
average of all aggregate sizes in the system. Both max- 
imal and average values increase constantly from 1.68 
(3.59) for 40/fc fc T to 18.04 (17.51) for 120fc 6 T. It becomes 
apparent, that for very weak hydrophobicities most of 
the surfactants are either isolated in bulk phase or in 
very small assemblies. For strong hydrophobicities, ag- 
gregates are very distinct and single surfactants in bulk 
phase are rare. There is, however, no sharp boundary or 
phase transition between small sub-micellar assemblies 
and proper micelles, but rather a continuous transition. 

Going back to the issue of nanocell division, we want 
to distinguish proper fission into nanocells of approxi- 
mately equal size from dissociation of sub-micellar ag- 
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FIG. 6 Mean aggregate and mean maximal aggregate num- 
bers for hydrophobicity values between 40 and 120fc(,T. The 
latter averages only the size of nanocells which are actually 
going to divide, while the former averages all aggregates. Sys- 
tem states have been averaged for 500r < t < lOOOr. Earlier 
states have been considered as transient. Below, parts of the 
final simulation states are shown for selected parameters. 



gregates. The previous analysis revealed that one cannot 
use absolute aggregate numbers, as they exhibit a trend 
for stronger hydrophobicity. Therefore, we characterized 
each event of the form of equation II II by the function 



Q(m, n) = 1 — 



\m — n\ 
m + n 



(16) 



which denotes the relative fission quality. Q(m,n) = 1 
for m = n, i. e. when the resulting nanocells are equal 
in size. Q becomes smaller as daughter cells become less 
alike. We have averaged Q as a function of the hydropho- 
bicity over all fission events in the simulation runs dis- 
cussed before. Mean and standard deviation are shown 
in figured (Q) varies between 0.41 and 0.65 with no sig- 
nificant trend for weak or strong hydrophobicities. More- 
over, standard deviation is very high. This reveals that 
fission into any two daughter cells is equally probable, no 
matter the ratio of their sizes. In terms of fission rates 
this finding can be written as 



k: 



k (m 



(17) 
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FIG. 7 Mean quality of nanocell divisions. A value of 1 
signifies division into daughter cells of equal size, while a value 
close to results from single surfactant dissociation. Error 
bars denote standard deviations. 
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FIG. 8 Number of surfactants as a function of time for dif- 
ferent catalytic rates k c . For slow metabolic turnover, the 
exponential shape of the auto-catalysis becomes apparent. 
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FIG. 9 Mean aggregate and mean maximal aggregate num- 
bers for catalytic rates between 0.03125 and 4. Or -1 . Consult 
the caption of figure|S|on how aggregation numbers have been 
obtained. 
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FIG. 10 Number of overall fission and fusion events (transi- 
tions of the form of equations 1131 and 1141 as a function of the 
catalytic rate k c . 



B. Influence of the catalytic rate on nanocell dynamics 

The catalytic rate has been varied from 2~ 6 t~ 1 to 
4.0t~ 1 in exponential steps. Hydrophobicity has been 
set to 80/cbT. Global surfactant dynamics are shown in 
figure|S| For slow metabolic turnover (< 0.125T" 1 ), the 
exponential shape of surfactant production becomes ap- 
parent, i. e. the constant precursor supply of precursors 
does not limit surfactant production over the simulated 
time span. This is tantamount to saying that unmetabo- 
lized precursor droplets are present throughout the whole 
simulation. Deceleration of the global dynamics is re- 
flected in slower fission rates of individual nanocells. For 
example, we observed that the first fission event is re- 
tarded by 180r on the average when catalytic rates are 
halved. For k c < 0.03125T" 1 , no fission occurs during 
the simulated time span. 

Figure shows mean aggregation numbers as a func- 
tion of the catalytic rate. As one can see, slow metabolic 
turnover increases both maximal and average aggrega- 
tion numbers (from 9.45 (5.27) for 0.5T" 1 to 18.55 (9.66) 
for 0.03125r _1 ). The trend becomes less apparent for 
faster metabolic turnover (1.0r _1 to 4.0V 1 ). For k c = 
0.125r , the maximal aggregation number is slightly 
higher than a monotonic trend would imply — a fact that 
we relate to statistical deviations, as fission events are 
considerably rare for small catalytic rates. The increase 
in both aggregation numbers is a natural consequence of 
the decelerated metabolism: when the precursor is me- 
tabolized slowly while its supply is held constant, the size 
of the hydrophobic core increases, and offers a bigger area 



for surfactants to attach. Hence, the maximal aggregate 
numbers increase. 

For the above runs, the number of fission and fusion 
events has been measured (see figure E3>- For all simula- 
tion runs, fusion and fission events are more or less bal- 
anced. This reveals that most of these events result from 
surfactant exchange with the bulk phase or from series of 
temporary fission and fusion during a single division pro- 
cess rather than from proper nanocell divisions. For low 
catalytic rates (0.03125r _1 < k c < 0.25r _1 ) the num- 
ber of such balanced transitions falls significantly from 
1028 for k c = 0.5-r- 1 to only 8 for k c = 0.03125t _1 . As 
in the case of aggregation numbers, the trend in transi- 
tion numbers can be related to the hydrophobic core: the 
more hydrophobic particles in the interior of a nanocell, 
the less dissociations occur on its surface. A strong anti- 
correlation between aggregation numbers and the num- 
ber of fission/fusion events (with a correlation coefficient 
of —0.917) justifies this hypothesis. 

It has to be pointed out, however, that the cat- 
alytic rate might affect nanocellular dynamics only dur- 
ing a certain transient time. It has been shown how 
the nanocellular dynamics depend on the ratio between 
metabolic turnover and precursor supply. Precursors are 
supplied by diffusion. Therefore, the rate of their incor- 
poration into an individual cell depends on the overall 
concentration of nanocells. Once a critical cell concentra- 
tion is reached, precursor incorporation might be slower 
than its metabolic turnover, which would undermine the 
above discussed effect. While the duration of such tran- 
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sient will depend on the ratio of precursor supply and 
turnover, dynamics after the transient might be little af- 
fected again. Ideally, simulations would be performed in 
a homeostatic system, with an influx of precursor solu- 
tion at one side and an outflux of reaction products at 
the opposite site of the system. Up to now, however, lit- 
tle is known about the performance of DPD in such open 
systems. 



IV. DISCUSSION 

In this paper we have presented an information-free 
nanocell based on a micellar system and a single auto- 
catalytic reaction that serves as metabolism. This simple 
system can be understood as a minimal self-replicating 
chemical system. As such, it denotes the boundary be- 
tween pure auto-catalysis and a more complicated self- 
reproducing system which would also include inheritable 
information. We analyzed the dynamics of this nanocell 
using a dissipativc particle dynamics approach. This sim- 
ulation technique can cover the relevant time scale, while 
it has been shown to be still physically accurate com- 
pared to other simulation techniques. As a consequence, 
we have been able to perform analyses of the system in 
a level of detail, that has — as far as we know — not been 
reached before in the study of self-replicating entities. 

The general replication cycle of micellar nanocells by 
metabolic turnover and division is very robust against 
changes in hydrophobicity and catalytic rates. It has 
been shown that the mean aggregation number of 
nanocells depends on the hydrophobicity of the surfac- 
tant (and precursor) as well as on the catalytic rate of the 
metabolism. For increasing hydrophobicity, a monotonic 
change in aggregation number with no sudden phase 
transition has been observed, ranging from a nearly ho- 
mogeneous solution with only submicellar aggregates for 
weak hydrophobicities to the formation of distinct mi- 
celles in surfactant-free water for a very high hydropho- 
bicity. The same monotonic increase in aggregation num- 
ber could be observed for increasing catalytic rates, i. e. 
fast metabolic turnover. 

It has been found that the rate of nanocell fission and 
surfactant dissociation depends on the size of the hy- 
drophobic core of the nanocells, and is more likely to oc- 
cur for small values in hydrophobicity and slow metabolic 
turnover. Daughter cells resulting from a fission event 
have been shown to vary significantly in size. There is 
neither a trend in the average size ratio of fission products 
nor in its variance. 

Our work shows that the envisioned replication-cycle 
of nanocells — namely incorporation and turnover of pre- 
cursor droplets followed by eventual aggregate division — 
is achievable over a wide range of parameters. In fact, 
there is no parameter combination for which the general 
replication cycle has been rendered impossible: although 
we have been able to decrease mean aggregation num- 
bers and increase dissociation rates up to a point where 



the system obeys no clear structures, we could not ob- 
serve that dissociation of single surfactants jeopardizes 
the growth of a otherwise stable nanocell. 

We have demonstrated the feasibility of a self- 
replicating system in the absence of genetic information. 
Although such a system is not able to evolve towards 
higher complexity, it could have served as a function- 
ing subsystem providing metabolism and embodiment for 
subsequent protocells of higher evolutionary complexity. 
It is conceivable, that independently evolved information 
systems like RNA might have become incorporated into 
such functioning replicators. When the two formerly in- 
dependent replication cycles of container and genome are 
orchestrated by coupling, such that each daughter cell of 
the dividing container is loaded with exactly one copy of 
the genomic information, one would obtain a true self- 
reproducing protocell with the ability to metabolize, di- 
vide and evolve. 

Apart from prebiotic scenarios in which micelles are 
considered as possible ancestors of subsequent vesicle 
based organisms, such systems are explored in the con- 
text of so-called living technology, i. e. artificial sys- 
tems that mimic life-like behavior. Our results suggest 
that the generic replication cycle of micellar nanocells 
is a robust basis for artificial life forms. We are cur- 
rently exploring one design of such an artificial protocell 
in which genomic information is coupled to a micellar 
self-replicating system similar to the one presented here. 
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